Silicon and III-V compound nanotubes: structural and electronic properties 
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Unusual physical properties of single-wall carbon nanotubes have started a search for similar 
tubular structures of other elements. In this paper, we present a theoretical analysis of single-wall 
nanotubes of silicon and group III-V compounds. Starting from precursor graphene-like structures 
we investigated the stability, energetics and electronic structure of zigzag and armchair tubes us- 
ing first-principles pseudopotential plane wave method and finite temperature ab-initio molecular 
dynamics calculations. We showed that (n,0) zigzag and (n,n) armchair nanotubes of silicon hav- 
ing n > 6 are stable but those with n < 6 can be stabilized by internal or external adsorption of 
transition metal elements. Some of these tubes have magnetic ground state leading to spintronic 
properties. We also examined the stability of nanotubes under radial and axial deformation. Owing 
to the weakness of radial restoring force, stable Si nanotubes are radially soft. Undeformed zigzag 
nanotubes are found to be metallic for 6 < n < 11 due to curvature effect; but a gap starts to open 
for n > 12. Furthermore, we identified stable tubular structures formed by stacking of Si polygons. 
We found A1P, GaAs, and GaN (8,0) single-wall nanotubes stable and semiconducting. Our results 
are compared with those of single-wall carbon nanotubes. 

PACS numbers: 73.22.-f, 68.43. Be, 73.20.Hb, 68.43.Fg, 61.46.+W 



-a 
c 

o 
o 



> 

(N 

in 
o 
in 
o 
-i— > 



i 

C 

o 
o 



X 



I. INTRODUCTION 

Carbon nanotubes are unique one dimensional 
nanostructurespj with their exceptional mechanical, elec- 
tronic and magnetic properties. 0, 01 01 01 While the 
use of single- wall carbon nanotubes (SWCNT) requires a 
completely new paradigms in the development of nanode- 
vices, Si still continues to attract interest for electronic 
applications in nanotechnology. Therefore, Si based nan- 
otubes have been subject of experimental and theoretical 
analysis. 

Even if a single-wall Si nanotube (SWSiNT) has never 
been observed, theoretical predictions have been per- 
formed for various kinds of Si tubes. Fagan et aL0, 
have investigated the structural and electronic properties 
of chiral SWSiNTs based on Density Functional The- 
ory (DFT). Barnard et ai.@| have examined the de- 
pendence of heat of formation and binding energy of 
SWSiNTs on their radius and chirality. The stability 
of (10,0) SWSiNT has been examined by using empiri- 
cal Monte-Carlo molecular dynamics method and found 
that it is stable at finite temperature. Q Ivanovskaya et 
aL0 investigated hypothetical Si nanotubes containing 
regular chains of metallocarbonedrenes using one dimen- 
sional tight-binding model within Huckel approximation. 
By using ab-initio calculations, Dumitrica et at |l0l| de- 
scribed how the smallest (2,2) and (3,0) SiNTs are sta- 
bilized by the axially placed metal atoms from differ- 
ent groups of periodic table. Ponomarenko et al.\n\ 
studied the energetics and relative stability of infinite 
and finite, clean and hydrogenated open-ended Si nan- 
otubes by using the extended Brenner potential. The 
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existence of H-doped stable tube-shaped finite SiNTs 
have been predicted 12] and their electronic structures 
have been compared with those of carbon nanotubes. 
Seifert et a/.jlj] have argued that structures of silicate 
and SiH nanotubes are more stable than bare Si nan- 
otubes. Singh et a/.[15| have investigated the stability 
of finite and infinite hexagonal prismatic structures of 
Si with 3d magnetic elements and predicted that such 
structures can be stabilized through doping by the ele- 
ments of transition metal (TM). Fullerene-structured Si 
tubulars, possibly based on Si24 have been produced. fl6| 
More recently, the successful synthesis of multiwalled Si 
nanotubes has been reported, [ljj Now, SWSiNTs are no 
longer hypothetical structures and it is not unrealistic 
to expect their fabrication with controllable size and di- 
ameter. Similarly, achievements of synthesis of single- 
wall BN nanotubes [3 and GaN [3123, A1N H3 thick- 
wall tubular forms has increased the interest in the the- 
oretical analysis of compound nanotubes. [22l l23l l2~iL l25| 
In addition, the synthesis of Mo and W chalcogenid 
nanotubes prj I2H |28| . and also NiCl tubular and cage 
structures have been realized. |29j 

In this paper we present a theoretical analysis of Si- 
and III-V compound nanotubes based on a state-of-the- 
art first-principle calculations. Our work is concentrated 
mainly on the tube structures which can be viewed as the 
rolling of graphene like honeycomb planes of Si or III-V 
elements on a cylinder of radius R. Starting from the 
precursor graphene-like honeycomb structures we investi- 
gated their stability, energetics and electronic properties 
of these nanotubes. Since O, O2, Si, Au and H are crit- 
ical elements for various processes on Si, we also exam- 
ined the adsorption of these atoms on SWSiNT. Finally, 
we studied the stabilization of unstable, small-diameter 
SWSiNTs through the internal and external adsorption 
of transition metal elements. In addition, we found that 
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tubular structures which are generated by stacking of 
triangles, pentagons and hexagons of Si are stable and 
metallic. The (8,0) zigzag tubes of A1P, GaN and GaAs 
are stable and semiconducting. The results obtained 
from the present study have been compared systemat- 
ically with those of SWCNT. The stable tube structure 
predicted in this study are hoped to motivate experi- 
mental research aiming at the synthesis of various tubu- 
lar structures of group IV elements and III-V and II- VI 
compounds 



II. METHOD 

We have performed first-pri ncip les plane wave 

calculations 30J |3l|| within DFT|32| using ultra-soft 
pseudopotentials. |3ll |3^| The exchange correlation po- 
tential has been approximated by Generalized Gradi- 
ent Approximation. (GGA) 34] Structures incorporat- 
ing TM atoms have been calculated using spin-polarized 
GGA. For partial occupancies we use Methfessel-Paxton 
smearing method. |35| The width of smearing is chosen 
between 0.01-0.1 eV depending on the system. All struc- 
tures have been treated by supercell geometry using the 
periodic boundary conditions. To prevent interactions 
between adjacent structures a large spacing (~ 10A) has 
been taken. Convergence with respect to the number 
of plane waves used in expanding Bloch functions and 
k-points in sampling the Brillouin zone are tested for 
the parent bulk crystals as well as tubular structures. 
In the self-consistent potential and total energy calcula- 
tions Brillouin zone of nanotubes has been sampled by 
(1x1x19 ) m esh points in k-space within Monkhorst-Pack 
scheme. |3qj Calculations of graphene and graphite struc- 
tures have been carried out using (19x19x1) and (8x8x6) 
k-point samplings, respectively. A plane-wave basis set 
with kinetic energy cutoff 200eV < h 2 \k + G\ 2 /2m < 
330eV has been used. All atomic positions and lattice 
parameters are optimized by using conjugate gradient 
method where total energy and atomic forces are min- 
imized. The convergence for energy is chosen as 10 -5 eV 
between two ionic step, and the maximum force allowed 
on each atoms is 0.05 eV / 'A. 

It should be noted that DFT based methods provide 
reliable predictions for the ground state properties , but 
band gaps are usually underestimated. Hence the energy 
band structure obtained from the single particle energy 
eigenvalues of Kohn-Sham equations are only approxima- 
tions to the real energy bands. Proper many-body self- 
energy corrections can be made by using GW method. |37| 
Recently, GW energies are compared with DFT-LDA 
results of (n, 0) SWCNTs which indicates shifts of va- 
lence and conduction bands and considerable increase of 
LDA band gap from 0.2 eV to 0.6 eV.[3£| It is suggested 
that GW corrections are small for large radius SWCNTs. 
Performing first-principles many-body Green's function 
calculations Spataru et aL[39j showed that the optical 
spectrum of both semiconducting and metallic small- 



radius SWCNTs exhibit important excitonic effects due 
to quasi-one dimensional nature. It is interesting to note 
that while the band gaps of (9,0), (12,0) and (15,0) zigzag 
SWCNTs have been measured by Scanning Tunneling 
Spectroscopy 01 to be 80, 42 and 29 meV, respectively, 
same band gaps have been predicted by GGA calculation 
piU to be 93, 78 and 28 meV, respectively. 

The stability of the structures we studied is the most 
crucial aspect of our work, since it provides valuable in- 
formation for the synthesis of these materials in future. 
In this respect an extensive analysis of stability has been 
carried out for various nanotubes. First, we applied a 
radial deformation to certain nanotubes and optimized 
their structures to see whether they relax to their origi- 
nal, undeformed circular forms under zero external force. 
Furthermore, we have performed, finite temperature ab- 
initio molecular dynamics calculations up to 1000 K us- 
ing Nose thermostatic for 250 time steps (0.5 ps) to 
check whether the optimized structure will be affected 
from random thermal motion of atoms or they maintain 
their tubular form at high-temperature. We believe that 
if there were any kind of structural instability it would 
be initiated and also enhanced within these time-steps at 
high temperature. 



III. HONEYCOMB STRUCTURE OF SILICON 
AND III-V COMPOUNDS 

One of the main difficulty for synthesizing Si nanotubes 
seems to be the absence of 2D silicon layer similar to the 
graphene structure of carbon. This is traced to the fact 
that in contrast to carbon, sp 3 -hybridization in Si is more 
stable than sp 2 -hvbridization. [43j In view of this situa- 
tion we examined whether the graphene-like 2D sheet of 
silicon can be stable. Two dimensional hexagonal lattice 
forming a honeycomb structure in the xy-plane has been 
periodically repeated along z-axis with 10 A spacing to 
minimize inter layer interactions. In order to reduce the 
effects of the constraints to be imposed by using the prim- 
itive unit cell we performed structure optimizations on 
the (2x2) cell in the xy-plane. Our calculations revealed 
that the planar structure (where all atoms lie in the same 
plane) is metastable, but it is buckled by 0.45 A relative 
vertical displacement of alternate atoms on the hexagons. 
The gain of energy upon buckling is 30 meV/atom. The 
binding energy is calculated to be 4.9 eV/atom which is 
0.6 eV lower than Si diamond structure and the average 
distance between nearest Si atoms is 2.2A. In the rest 
of the paper this graphene-like structure will be speci- 
fied as buckled honeycomb structure. As shown in Fig^ 
the detailed band structure and total density of states 
(TDOS) analysis indicate that both buckled and planar 
system have large band gaps along TK and MK direc- 
tions, but conduction and valence bands cross the Fermi 
level at the k-point of Brillouin zone. The electronic 
structure of the system does not change significantly as 
a result of buckling, except some of the bands split due 
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Planar Buckled 




FIG. 1: Band structure and total density of states (TDOS) 
analysis for planar and buckled honeycomb structure (HC) of 
Si. Light-lines correspond to planar (having P6/MMM sym- 
metry) and dark-lines correspond to buckled structure (hav- 
ing P-3M1 symmetry). The zero of energy is set to the Fermi 
level Ef- The inset shows the 1/12 of hexagonal Brillouin 
zone. 



to the lowering of the rotation symmetry. Using a simi- 
lar method but different pseudopotentials and exchange 
correlation potential Takeda et al. |44| have examined pla- 
nar and buckled honeycomb structures of Si. Our results 
obtained in four times larger cell hence allowing more 
variational freedom are in overall agrement with the re- 
sults in Ref. |44|. Moreover, we performed an ab-initio 
molecular dynamics calculations on 2x2 supercell provid- 
ing further evidence that buckled honeycomb structure is 
stable at 500K for 250 time steps. 

Next we address whether a graphite like structure of 
Si (or graphitic Si) can form? Our study distinguished 
chemisorption and physisorption states in the interlayer 
interaction, in contrast to only the physisorption state 
in graphite. 45] The chemisorption state corresponding 
to a smaller lattice parameter c = 6A is energetically 
more favorable, namely the binding is 5.1 eV which is 
0.4 eV smaller than that of the bulk Si. We note, 
however, that the lattice parameters, binding energies 
depend on the approximation of exchange-correlation 
potential, [i^. |47| Present GGA calculations yield rela- 
tively larger interlayer distance and requires the incorpo- 
ration of long range Van der Walls (VdW) interaction. 
Such an analysis related with the binding of graphite 
showing that the lattice parameter calculated by GGA is 
improved upon including VdW attraction has been pre- 
sented elsewhere. |48^ 

Similar to Si, the honeycomb structures of A1P, GaAs 
and GaN are found to be also stable, but less energetic 
relative to the bulk crystal by 0.8, 1.1 and 0.6 eV per 
basis, respectively. However, the buckling is not favored 
in order to hinder the formation of dipole layer. 



IV. SINGLE WALL SILICON NANOTUBES 
A. Energetics and Stability 

Having discussed the stability of Si buckled honeycomb 
structure (Si-HC), now we present our systematic anal- 
ysis of (n, 0) zigzag and (n, n) armchair SWSiNTs for 
different n values; namely n = 3 — 14 for zigzag and 
n = 3,6,9 for armchair structures. The (3,0) zigzag 
SWSiNT has clustered upon structure relaxation indi- 
cating that it is not stable even at T=0 K. While the 
structure optimization has resulted in a regular (4,0) and 
(5,0) tubular structures, the ab — initio MD calculations 
showed that these nanotubes eventually transform into 
cluster at higher temperatures as shown in Fig. [2J Sig- 
nificant distortions can be easily noticed in (6,0) and 
also (7,0) SWSiNTs, but tubular character and hexag- 
onal structures on the surface have remained. The (6,0) 
zigzag tube, which has a radius of R = 3.8A as well as 
those with larger radii remain stable at temperatures up 
to 800 K. Barnard et al. also reported the instability 
of (3,0) SWSiNT in their first-principles study, but they 
considered (4,0) and (5,0) SWSiNTs as stable structures 
depending on their geometry optimization performed at 
T=0 K. Present results set a limit for fabricating small 
radius SWSiNTs. The first and second nearest neighbor 
interactions between Si atoms become relevant for the 
stability of small radius nanotubes and causes cluster- 
ings, if R < 3.8A Similar behavior is also obtained for 
(n, n) armchair SWSiNTs. For example (3,3) SWSiNT 
is clustered at 800 K in spite of the fact that geometry 
optimization yields tubular structure at T=0 K. On the 
other hand, the (6,6) tube having relatively larger radius 
remained stable at 800 K after 250 time steps. In contrast 
to (n, n) SWSiNTs which are found unstable for n < 6, 
the (3,3) SWCNT is known to be stable and experimen- 
tally fabricated. ^3113 The difference in the chemical be- 
havior of C and Si can be traced to the difference in their 
7r-bonding capabilities. Si tends to utilize all of its three 
valence p-orbitals, resulting in sp 3 -hybridization. In con- 
trast, the relatively large promotion energy from C-2s to 
C-2p orbitals explains how carbon will activate one va- 
lence p-orbital at a time leading, in turn, to sp, sp 2 , sp 3 - 
hybridizations in ID, 2D and 3D structures. This is the 
explanation why tubular structures of C are more stable 
than those of Si. Moreover, since the interatomic dis- 
tance increases significantly in going from C to Si, the 
7r — 7r overlap decreases accordingly, resulting in much 
weaker 7r-bonding for Si tubes in comparison with that 
for carbon tubes. 

After the discussion of stability, we next analyze the 
energetics, namely the behavior of binding energy (Ei,) 
as a function of the radius (or n) of the tube. per 
atom is calculated using the expression, 

E b = {Et[SW SiNT] - N(E T [Si])}/N (1) 

in terms of the total energy of the optimized SWSiNT 
having N Si atoms per unit cell, Et [SWSiNT] , and the 
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FIG. 3: The calculated binding energy per atom for (n, 0) 
zigzag and (n, n) armchair SWSiNTs. Dark region indicates 
that tubular structures are unstable at finite temperature. 
The calculated strain energies, Es per atom for (n, 0) and 
(n, n) SWSiNTs are shown by inset. 




FIG. 2: Structures of (4,0), (8,0) and (6,6) SWSiNTs at T=0 
and T=500K after 250 time steps. Tubular structure has 
remained in (6,6) and (8,0) SWSiNT, but (4,0) structure has 
clustered. 



Barnard et al. Calculated value of the strain energy of 
a zigzag SWSiNT is smaller than the strain energy of a 
zigzag SWCNT having comparable radius. 41] In classical 
theory of elasticity the strain (or curvature) energy of a 
tubular structure is given by the expression E$ — a/ R 2 , 
where a is a function of Young's modulus and thickness 
of the tube's wall.0, |H2 Result of the present calcula- 
tions in Fig|3gives a fair fit to the expression, a/R 2 with 
a ~ 2.07eV/A 2 . 



total energy of N, free Si atom Et[Si\. It is found that 
E\j ~ 4.9 eV and slightly increases as the radius R (or 
n) increases for both zigzag and armchair SWSiNTs as 
displayed in Fig|3J The energy increase with n is small. 
According to our results E^s of (n, n) armchair SWS- 
iNTs are ~ 0.05 eV larger than those of (n,0) zigzag 
ones because of their relatively larger radius at a given 
n. Corresponding E b for SWCNTs is calculated to be 9.1 
eV|5l| theoretically. 

Finally, the strain energy per atom is calculated rela- 
tive to the energy of honeycomb structure, 

E s = E b [SWSiNT] - (E b [Si - HC] ) (2) 

by subtracting the binding energy (per atom) of opti- 
mized honeycomb structure, E b [Si — HC] from the bind- 
ing energy of SWSiNT. Slight increase in strain energy 
is observed as the radius R or n decreases. This is 
an expected result, since the structure becomes more 
graphene-like with the increasing radius. Calculated 
strain energies given by inset in Fig|3| are also in agree- 
ment with the results obtained by Fagan et a/.Q and 



B. Mechanical Properties 

Radial flexibility is a criterion for the stability of tubu- 
lar structure. SWCNTs are known to be flexible for the 
deformations in radial directions; they can sus- 

tain severe radial deformation transforming the circular 
cross section into an elliptical one with minor (b) and 
major (a) axis. The radial deformation on a bare tube 
of radius R is specified in terms of the strain associ- 
ated with the pressing of the tube along the minor axis, 
e yy = (b — R)/R and the strain associated with the ex- 
pansion of the tube along the perpendicular major axis 
txx — ( a ~ R)/ R- Theoretical and experimental research 
have shown that radially deformed tubes relax reversibly 
to original circular cross section whenever the external ra- 
dial force is lifted. Moreover, radial deformation can 
modify the electronic structure reversibly, that leads to 
a tunable band gap engineering. [55l |5|J For example, a 
semiconducting (n,0) can be metallic under radial defor- 
mation. Our results indicate that SWSiNTs display a 
behavior different than that of SWCNTs. We performed 
a systematic analysis of radial strain for (8,0) zigzag and 



(6,6) armchair SWSiNTs. First, these tubes have been 
deformed by applying e yy —-0.1 r 0.2 and -0.3. Then the 
stress (or constraint) imposing these radial strains has 
been lifted and the structure has been optimized. Con- 
trary to situation in carbon nanotubes, up to the ap- 
plied strain e yy <-0.2 the SWSiNTs have remained in 
deformed state. For example, (8,0) tubes with initial ra- 
dial strain of e yy =-0.1 and -0.2 are relaxed to a plastic 
deformation corresponding to e yy =-0.09 and e OT =-0.14, 
respectively. Similar results have been obtained for (6,6) 
armchair SWSiNT with initial radial strain of e yy =-0.l 
and -0.2. In contrast, the tubes, which initially strained 
by e TO =-0.25 and -0.3 have relaxed to a state with negli- 
gible residual strain. The total energy of the undeformed 
SWSiNT have been found to be lower (more ener- 
getic) than the total energy E^(e yy ) of tubes which were 
relaxed upon radial deformation —0.3 < e yy < 0. How- 
ever, the energy difference AE = EJp(e yy ) — E^ > is 
very small. The weakness of 7r-bonds of Si as compared to 
carbon nanotubes is possibly a reason why the restoring 
forces are not strong enough to derive the deformed state 
to relax back to the original undeformed state. Once the 
applied radial deformation gets significant (\e yy \ > 0.2) 
the restoring forces become strong enough to derive the 
relaxation towards circular cross section. On the other 
hand, after a severe radial strain that causes to a signifi- 
cant coupling between opposite internal surfaces the de- 
formed state may be more energetic (i.e Ej.(e yy ) < E^ or 
it may relax to different structures such as clusters. This 
situation constitutes an important difference between Si 
and C single-wall nanotubes. 

Axial strength of SWSiNT, or the elastic stiffness along 
tube axis is defined as the second derivative of the strain 
energy per atom with respect to the axial strain e zz , 
namely k = d 2 Et / l de 2 zz . The elastic stiffness of the (8,0) 
SWSiNT along its axis is calculated to be 23 eV. This 
value is significant, but smaller than that of SWCNT 
which is calculated to be 52-60 eV.|57j 



C. Electronic Structure 

A systematic analysis of electronic structure indicates 
that metallic zigzag SWSiNTs 6 < n < 11 have three 
bands crossing the Fermi level, but a band gap between 
valence and conduction bands opens when n > 12. Sim- 
ilar effect has been obtained for zigzag SWCNTs when 
n > 7.0,^3 This metal-semiconductor transition was at- 
tributed to the energy shift of the singlet 7r*-band which 
is normally empty, but becomes filled due to increased 
fj* — 7r* hybridization at small radius. 0, In the 
present case it appears that a* — tt* hybridization be- 
comes significant at relatively larger radius. The con- 
ductance of all these infinite, perfect tubes (6 < n < 11) 
is predicted to be equal to 3G (G a = 2e 2 /h). Simi- 
lar metallic behavior is also obtained for armchair types 
namely for (6,6) and (9,9) SWSiNTs. The conductance 
of ideal infinite (n, n) tubes is 2G„, but not 3G as in 



5 



W 7 

8& / (7.0) 




(10,0) 












^^^^^^^ 


r z 


r z 


r z 



FIG. 4: The energy bands calculated for (7,0), (9,0), (10,0), 
(11,0), (12,0) and (14,0) SWSiNTs using GGA. The lowest 
conduction band, or singlet 7r*-band is indicated. The zeros 
of energy are set at the Fermi level. 

metallic (n,0) zigzag tubes. Fig. 0] presents the system- 
atic analysis of (n, 0) tubes for 7 < n < 14 and clearly 
shows how the singlet 7r*-band gradually raises as R in- 
creases. 

Based on LDA calculation Fagan et aL^,0| also found 
(6,6) and (6,0) SWSiNT's metallic, but they predicted 
(10,0) and (12,0) zigzag nanotubes are semiconductor 
with a small band gap of 0.1 eV. The disagreement be- 
tween the present one and those of Fagan et ai.0, 
may be due to the differences in pseudopotentials and 
in the approximation of exchange correlation potential. 
Note that the transition of (n, 0) SWSiNTs from metal- 
lic to semiconducting state through gap opening may oc- 
cur at n that is smaller than predicted by the present 
study as well by Fagan et ai, if self-en ergy correction 
are taken into account by GW method. |37| Electronic 
structure analysis performed for the tubes under strain 
both radially and axially showed that metallic character 
is not altered but only the position of Fermi level slightly 
changed due to deformation. The modification of elec- 
tronic structure with chirality may offer the possibility 
of fabrication of nanodevices using SWSiNT junctions. 
On the other hand, SWSiNTs can be used as metallic 
interconnects, since their conductance is not severely af- 
fected by deformation. 



D. Interaction of SWSiNT with Atoms and 
Molecules 

The interaction of Si nanotubes with oxygen atom and 
oxygen molecule is extremely important for technologi- 
cal applications. The adsorption of oxygen atom is stud- 
ied by placing it initially above a Si-Si bond parallel to 
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FIG. 5: (a) The optimized geometry of O atom adsorbed 
over the axial site of (8,0) SWSiNT. (b) Initial adsorption 
geometry of O2 over the axial site, (c) Dissociation of O2 
after relaxation and formation of Si-O-Si bridge bands over 
the adjacent zigzag Si-Si bonds, (d) Variation of interaction 
energy as a function of the distance d of O2 molecule from the 
parallel axial Si-Si band, (e) Optimized geometry of a single 
Si atom adsorbed on the top site (indicated by a dark ball). 



the axis of a (8,0) SWSiNTs. The optimized structure 
shown in Fig. [SJi. has a very strong chemical bonding 
between O and SWSiNT with E b =8.1 eV and the nearest 
Si-0 distance 1.7 A. The resulting geometry showed that 
SWSiNT is slightly distorted upon O adsorption. 

The interaction between O2 and SWSiNT has been re- 
vealed by calculating the binding energy as a function of 
the separation d from the axial Si-Si bond of the tube as 
shown in Fig. |SJb). O2 molecule is kept unrelaxed and 
taken parallel to the Si-Si bond. The calculated energy 
versus distance curve E(d) in Fig. Eld) shows that O2 
can be attracted to the tube, but there is no physisorp- 
tion state as in 2 +SWCNT.[48| The minimum of E{d) 
occurs at 1.9 A. Upon relaxation of the tube and O2 near 
this minimum, the molecule has been dissociated to form 
two Si-O-Si bridge bonds over the zigzag Si-Si bonds, 
and concomitantly SWSiNT has been distorted locally 



as illustrated in Etc). The distances between nearest Si- 
O and 0-0 are 1.7 and 3.3 A, respectively. We repeated 
the structure relaxation by initially placing O2 at a larger 
distance d — 2.hA from the surface of the tube and we 
obtained the same dissociated state. Our results indi- 
cate that there will be a strong interaction between Si 
nanotube and oxygen molecule in open air applications. 

The SWSiNT surface is found to be reactive against 
Si, H and Au atoms. Si atom attached to the top site 
is bound by Ef, ~ 5 eV. One Si atom of the tube is 
plunged inside the tube and a small cluster is formed 
at the surface (see Fig. OJe)). The chemisorption energy 
of H and Au atoms is strong and found to be 4.4 eV and 
3.4 eV, respectively. 



V. STABILIZATION OF SILICON NANOTUBES 
BY TRANSITION METAL ATOM DOPING 

Recently, Singh et al. 15] showed that Si-clusters and 
Si-tubular structures formed by top-to-top stacking of Si- 
hexagons can be stabilized by the implementation of TM 
atoms inside these structures. Those structures are not 
only stabilized, but also acquired magnetic properties. 
Earlier, TM atoms are shown to form rather strong bonds 
with the carbon atoms on the surface of SWCNTs.jSlJ 
Motivated by the work of Singh et a/.0], we investi- 
gated whether (3,0) SWSiNT can be stabilized in the 
same manner. The (3,0) tube has radius R ~ 2AA in 
which interaction between the atoms located at the op- 
posite walls of the tube as well as excess strain on the 
Si-Si bonds are the prime causes of structural instability 
even at T=0 K. On the other hand, the radius of (3,0) 
is comparable with the sum of ionic radii of V and Si, 
i.e Ry + Rsi — 2.27 'A, and hence V atoms can easily be 
accommodated inside the tube. We considered a (3,0) 
SWSiNT, which has V atoms implemented inside and 
periodically arranged along the tube axis. Because of su- 
percell geometry used in the calculations both chains of 
V atoms (V-LC) and (3,0) tube have common lattice pa- 
rameter (See Fig|B}. The optimized structure, that con- 
sists of planar hexagons are stacked with V-LC passing 
through their centers, has been found to be stable. The 
energy of the V-stabilized structures is lowered by 12.9 
eV relative to the energies of the V-LC and Si-tube with- 
out V-LC in it but having the same atomic structures 
as Si tube with V-LC. Spin-relaxed calculations resulted 
in zero magnetic moment fj, = 0. The p-d hybridization 
between Si and V atomic orbitals is the cause of stability 
and lowering the total energy. 

The radius of the unrelaxed (3,3) SWSiNT of 3.7 A 
is too large and may not be suitable for its stabilization 
through the implementation of an atomic chain. Per- 
haps, it may be better suited to accommodate a small 
clusters of atoms. We considered the possibility whether 
the (3,3) SWSiNT is stabilized by TM atoms adsorbed 
on the external surface of the tube. To this end, we 
studied Ti and Cr atoms adsorbed on the hollow sites 



(3,0) + V 



(3,3) + 3Cr 




(3,3) + 3Ti 



(3,3) + 6 Cr 



FIG. 6: Perspective and top (cross section) view of the fully 
relaxed V-doped (3,0), Ti and Cr doped (3,3) SWSiNTs. Our 
results revealed that small radius SWSiNTs can be stabilized 
by doping transition metal element. 
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(z.e above the hexagons formed by Si atoms). As shown 
in Fig. we examined 3 different (3,3)+TM structures; 
namely 3 Ti-LC, 3 Cr-LC and 6 Cr-LC are adsorbed on 
the (3,3) SWSiNT surface. Spin-relaxed GGA calcula- 
tion are carried out to optimize the geometric structure. 
The external absorption 3 Ti-LC (or 3 Ti atom per unit 
cell of (3,3) tube) prevented the tube from collapsing 
into a cluster, but the circular cross section changed to 
a polygonal one. The ground state has been predicted 
to be non-magnetic with /i = 0. The external adsorp- 
tion of 3 Cr-LC also resulted in polygonal cross section, 
but ferromagnetic ground state with net magnetic mo- 
ment = 9.7 [Ib (Bohr magneton). The circular cross 
section is maintained by the adsorption of 6 Cr-LC. This 
latter structure has also ferromagnetic ground state with 

n = yi.2hb- 

Calculated energy band structure of (3,0)+V, 
(3,3)+Ti, (3,3)+3Cr and (3,3)+6Cr are presented in Fig. 
13 The (3,0)+V structure is a metal. Six bands crossing 
the Fermi level yields quantum ballistic conductance of 
G=6Go- The partial density of states indicates that V- 
3d and Si- 3p orbital character dominate the states at the 
Fermi level. The (3,3)+3Ti structure is a semiconductor 
with a very narrow band gap. In the case of (3,3)+3Cr 
and (3,3)+6Cr several majority (spin-up f) and minority 
(spin-down j) bands are crossing the Fermi level. Hence 
both structures are metals with finite density of majority 
D(E = Ep, t) and minority D(E = Ep, {) spin states at 
E F . However, D(E = E F , |) - D(E = E F , j) is signif- 
icant. These properties, which are also depends on the 
decoration of the tubes, can be used in nanospintronic 
device applications. Much recently Dumitrica et aL|lO| 



FIG. 7: The energy band structures of (3,0)+V, (3,3)+3Ti, 
(3,3)+3Cr and (3,3)+6Cr calculated by using SCF spin po- 
larized GGA method. Solid lines and dotted lines are are for 
majority and minority states, respectively. The inset shows 
the density of majority (dark) and minority (light) spin states 
at the Fermi level of (3,3)+3Cr and (3,3)+6 Cr structure. 



have investigated the stabilization (3,0) zigzag and (2,2) 
armchair SWSiNTs by various atoms (Zr, Sc, Ti, Cr, Fe, 
Ni, Be, And Co) axially placed inside the tube. However, 
they consider neither the magnetic ground state due to 
specific TM atoms, nor the stabilization of tubes having 
relatively larger radius. 



VI. OTHER TUBULAR STRUCTURES OF 
SILICON 

We now consider different tubular structures which are 
formed by stacking of atomic polygons as presented in 
Fig. [S] Three tubular structures given in the first row, 
i.e. Ti, T2 T3 are made of triangles of silicon atoms, 
which are stacked along the axis of the tube. In Ti , tri- 
angles are identical and placed in top-to-top (or eclipsed) 
position; in T2 the triangles are staggered; T3 structure 
is constructed by insertion of a Si-LC into the Ti struc- 
ture as such that the chain atoms are centered in between 
the layers. The same convention is followed in labelling 
the tubular silicon nanowires with pentagonal (Pi, P2, 
P3) and hexagonal (Hi, H2, H3) cross-sections. 

Starting with those geometries, structural optimiza- 
tions in large supercells yielded T2, T3, and H3 as unsta- 
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TABLE I: Structural and conductance properties of silicon 
tubular structures that are found stable. Et is the binding 
energy per atom, c is the unit cell length of the periodic struc- 
ture, d is the in-plane bond lengths of the polygonal atomic 
layers. Equilibrium conductance values are given by G in 
units of conductance quantum, G = 2e 2 jh 



Structure 


E b (eV) 


c(A) 


d(A) 


G (2e 


Ti 


4.62 


2.37 


2.38 


6 


T 2 










T 3 










Pi 


4.79 


2.37 


2.37 


10 


P 2 


4.68 


2.70 


2.46 


10 


Pa 


4.66 


4.10 


2.58 


6 


Hi 


4.77 


2.29 


2.37 


6 


H, 


4.74 


2.52 


2.42 


9 


H 3 











FIG. 8: Some tubular structures of Si. The tubes are formed 
by stacking of Si polygons; the labels T, P, and H stand for 
triangular, pentagonal, and hexagonal configurations, respec- 
tively. Ti, Pi, and Hi tubes have top-to-top arrangement of 
layers. In T2, P2, and H2 the layers are staggered. T3, P3, 
and H3 structures have extra atoms centered in between the 
layers. 




FIG. 9: Energy band structure of Si tubular structures formed 
by top-to-top stacking of triangular (Ti), pentagonal (Pi), 
and hexagonal (Hi) polygons. In the right panels the corre- 
sponding equilibrium conductance in G = 2e 2 /h curves are 
plotted. The Fermi levels are set to zero. 



ble, which deformed into clusters. The structural param- 
eters, the binding energies, and equilibrium conductance 
values of stable tubular structures are summarized in Ta- 
ble [I] All the stable silicon tubes have comparable bind- 
ing energy values, with Pi structure being slightly more 
favorable energetically. Also within each set of tubular 
structures, the top-to-top arrangement provides the high- 
est cohesive energy. In Fig. [5] the energy band structure 
and equilibrium conductance plots of T l5 Pi, and Hi 
structures are presented. All of them are found metallic 
with calculated ideal conductance values of 6G0, 10Go, 
and 6G0, respectively. A common feature noticed in the 
energy band structures of Ti, Pi, and Hi tubes is that 
there are almost filled bands available in close vicinity of 
the Fermi levels which may lead to drastic conductance 
variations due to small structural perturbations in the 
tubes or small bias voltage. 

Earlier, the pentagonal nanowires of Si and also those 
of several metals such as Na, Al, Cu, Pb, Au, Fe, Ni, 
and Xe were investigated by Sen et al. j5^ using similar 
calculation methods. The results for a specific structure, 
P3, which is common in both studies, is in agreement. 
The stability of the P3 is further strengthened by the 
finite-temperature ab-initio molecular dynamics calcula- 
tion carried out in the present study. Recently, Bai et 
al. [60| have studied the stability of infinite and finite Si 
(top-to-top square), Pi and Hi structures by performing 
ab-initio calculations using different methods including 
pseudopotential plane-waves and classical molecular dy- 
namics calculations at 2000K. Their stability analysis for 
Pi and Hi are in agrement with present results. In addi- 
tion to Pi and Hi, the present study deals with Ti and 
also staggered ones. 
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VII. SINGLE WALL NANOTUBES OF III-V 
COMPOUNDS 

Motivated by interesting properties of SW(BN)NT and 
opto-electronic and field emitting properties of GaN and 
A1N tubular forms [H HI we choose (8,0) AlP, GaAs 
and GaN single-wall nanotubes as prototype to investi- 
gate the stability and electronic properties of III-V com- 
pound nanotubes. Even if the single-wall nanotubes of 
these compounds have not been sythesized yet, the pre- 
dictions of present work is essential for further efforts to 
achieve it. The initial bond lengths are chosen as the 
distance between nearest cation and anion atoms in bulk 
structure. After relaxation of all atomic positions, as 
well as lattice constant c, the tubular structures remained 
stable. The ab-initio MD calculations also showed that 
SW(A1P)NT remained stable at room temperature after 
250 time steps. Eb is calculated to be 9.6 eV per AlP 
basis. The radius of the tube is 5.2 A. The structure 
is not a perfect tube but the hexagons on the surface 
are buckled. The nearest Al-P distance is 2.3A, and sec- 
ond nearest neighbor distance i.e, nearest P-P and Al-Al 
distances are 3.9A and 3.8A, respectively. The energy 
band and TDOS analysis in Fig. ^Jpoints out that (8,0) 
SW(A1P)NT is a semiconductor (insulator) with a band 
gap of 2.0 eV. 

Initial tubular structure of (8,0) SW(GaAs)NT is 
maintained after geometry optimization at T=0 K. Sim- 
ilar to SW(A1P)NT, hexagons are buckled. Eb is calcu- 
lated to be 7.7 eV per GaAs and the radius is 4.8 A. The 
nearest Ga-As distance is 2.4A, and nearest Ga-Ga and 
As- As distances are 3. 9 A and 4.1 A, respectively. The 
(8,0) SW(GaAs)NT is also a semiconductor (insulator) 
with a band gap of 0.9 eV. 

We place a special emphasis on GaN nanotubes. 
which are successfully synthesized by an epitaxial casting 
method where ZnO nanowires are initially used as tem- 
plates. GaN nanotubes produced this way have diameter 
of 300 A and minimum wall thickness of 50 A. They 
are semiconducting and hence they would be a possi- 
ble candidate for optoelectronic applications. Whether 
single-wall GaN tube of smaller diameter (2R ~ lOA) 
can be stable and can exhibit technologically interesting 
electronic properties is important to know. We again 
took (8,0) SW(GaN)NT as prototype for the sake of con- 
sistency. Stable tubular geometry is obtained by both 
geometry optimization at T=0 K and ab-initio MD anal- 
ysis at T=800 K. Upon relaxation atoms on the sur- 
face are buckled. Eb is calculated to be 11.5 eV per 
GaN and the radius is 4.1 A. The nearest Ga-N dis- 
tance is 1.8A, and nearest Ga-Ga and N-N distances are 
3.1A and 3.2A, respectively. We found that the (8,0) 
SW(GaN)NT is a semiconductor (insulator) with a band 
gap of 2.2 eV. Previous first-principles study performed 
by Lee et al.[2j| by using LDA method predicted the 
similar band gap for SW(GaN)NT. Here, we examine 
also whether SW(GaN)NT is radially elastic. To this 
end we started with the elliptically deformed nanotube 
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FIG. 10: Energy band structures (left panels) , total density of 
states (TDOS) and partial density of states (PDOS) on atoms 
(right panels) of (8,0) zigzag SW(A1P)NT, SW(GaAs)NT, 
and SW(GaN)NT. Anion (Al, Ga) and cation (P, As, N) con- 
tributions to TDOS are shown by dashed and light-continuous 
lines. Zero of the energy is set to the Fermi energy. 



under e yy — —0.1, and let it relax in the absence of radial 
forces. Similar to the Si nanotube, SW(GaN)NT is found 
to be radially soft 

For the sake of comparison, the binding energy and 
bond distance of bulk crystal, honeycomb structure, and 
(8,0) tubular structures of Si, AlP, GaN, and GaAs are 
presented in Table ITT1 In these structures the covalent or 
(covalent-ionic) mixed bonds have different orbital com- 
binations. While bulk crystals are tetrahedrally coordi- 
nated and have bonds formed by sp 3 -hybrid orbitals, in 
honeycomb and tubular structures bonding through sp 2 - 
hybrid orbitals dominates the cohesion. 
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Si 


A1P 


GaN 


GaAs 




B H T 


B H T 


B H T 


B H T 


E„ (eV) 

d(A) 


5.4 4.9 4.8 
2.3 2.3 2.3 


10.4 9.6 9.6 
2.4 2.2 2.3 


12.4 11.8 11.5 
2.0 1.8 1.8 


8.3 7.2 7.7 

2.4 2.2 2.4 



TABLE II: Calculated binding energies Et and bond distances 
d in A of various structures of Si, A1P, GaAs and GaN. B:3D 
bulk crystal, H:2D honeycomb structure, T:(8,0) single wall 
nanotube structure. The units of binding energy is eV/atom 
for Si, and eV/basis for III-V compounds. 



VIII. CONCLUSION 

In this paper, we analyzed the stability of Si and III- 
V compound, single-wall nanotubes and calculated their 
optimized atomic structure and energy band structure. 
Si as well as III-V compounds can form stable 2D hon- 
eycomb structure, which is precursor of nanotubes. The 
energy necessary to roll these honeycomb structures over 
a cylinder of radius R to make a perfect nanotube is how- 
ever small as compared to those in carbon nanotubes. 
We found that Si single-wall nanotubes with small ra- 
dius are unstable and are clustered either at T=0 K or 
at finite temperatures. For example, while (3,0) is un- 
stable even at T=0 K, (4,0) and (5,0) lose their tubular 
character and tend to form cluster at T=500 K. Stable 
(n, 0) zigzag SWSiNTs are metallic for 6 < n < 11, but 
become semiconducting for n < 12. The metallicity of 
small radius (n, 0) tubes is a typical curvature effect and 
is resulted from the dipping of the singlet 7r*-band into 
the valence band at small radius. Stable (n,n) armchair 
SWSiNTs (n=6,9)are metallic. Our study on radially 
deformed (8,0) and (6,6) SWSiNTs demonstrated that 
these nanotubes are radially "soft", and hence are de- 



void of strong restoring force that maintains radial elas- 
ticity. The radial softness of Si tubes is a behavior which 
distinguishes them from carbon nanotubes. In contrast 
to that axial stiffness the Si nanotube has ben found to 
be high. We predicted that oxygen molecule adsorbed on 
the Si-Si bonds dissociates. A strong interaction between 
O/O2 and SWSiNT appears to be serious in future pro- 
cesses involving Si tubes. Adatoms like Si, Au and H can 
also form strong chemisorption bonds with the atoms on 
the surface of SWSiNT. We showed that unstable, small 
radius SWSiNTs can be stabilized through the imple- 
mentation or external adsorption of 3d transition metal 
atoms. In particular, the decoration of the tube surface 
by the external adsorption of transition metal atoms can 
lead magnetic properties which may find potential tech- 
nological applications. Small radius tubular structures 
different than those based on honeycomb structure have 
been identified. Finally, we found III-V compound (8,0) 
nanotubes (A1P, GaAs, and GaN) stable at least at room 
temperature and they are semiconductor with band gap 
ranging from 0.9 eV to 2.2 eV. In contrast to small ra- 
dius metallic Si nanotubes, (8,0) compound nanotubes 
are semiconductor. The band gap increases with de- 
creasing row number of elements. Even though not all 
the structures treated in this study have not been real- 
ized experimentally yet, the predictions obtained from 
the present first-principles calculations are expected to 
be essential for further research in this field. 
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